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Abstract 

The production of the heavy charmoniuni and (c6)-quarkonium, through the top quark semi- 
exclusive decays via the flavor changing neutral currents (FCNC), has been systematically studied 
within the non-relativistic QCD. In different to the conventional squared amplitude approach, 



Ph. 

Q^. to simplify the results as much as possible, we adopt the 'improved trace technology' to do our 

(D ■ 
j^, '. calculation, which deals with the hard scattering directly at the amplitude level. If assuming 

the higher excited heavy-quarkonium states, such as the color-singlet and spin-triplet S'-wave 

> : 

C^ ■ state Ip^iji), the color-singlet P-wave states |[^Pi]i) and ipPj]!) (with J = 0,1,2), and the two 

O , 

cn , color-octet components |[^5o]85') and |['^S'i]8(7), decay to the ground color-singlet and spin-singlet 

■r^lj- I S-wave state |[^5'o]i) with 100% efficiency via the electromagnetic or hadronic radiations, we 

cn ■ obtain the total decay width: Tt_,\(^^^)[iSo]^) = ITl.l+gg" KeV and ^t~,\(cb)[^So]i) = '^■^'^^i.il KeV, 

where the uncertainties are caused by varying m^ = 172.0 it 4.0 GeV, rrif, = 4.90 it 0.40 GeV and 



^ I rric = 1.50 lb 0.25 GeV. At the LHC with the center-of-mass energy \/S = 14 TeV and a high 

luminosity C oc 10^^cm~^s~^, sizable heavy-quarkonium events can be produced through the top 
quark decays via FCNC. 
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I. INTRODUCTION 

With a high luminosity C oc 10^^cin~^s~^ and a high proton-proton coUision energy 
\fS = 14 TeV, a large amount of t-quark events about the order of 10^ per year will be 
generated at the large hadronic collider (LHC) [ll, l2l- This makes LHC a 'top factory' that 
shall provide us a unique platform to study the heavy quarkonium properties in comparison 
to the proton-anti-proton collider TEVATRON and the e~^e~ colliders such as Belle, BABAR, 
BEPC, CLEO-c and etc.. 

Being the heaviest fermion with a mass close to the electroweak symmetry breaking scale 
in the standard model (SM), the top quark is speculated to be a sensitive probe of new physics 
beyond the SM. In the literature, many studies have shown that the top quark rare decays 
via a flavor changing neutral current (FCNC), i.e. t — )■ Vq, {V = 'y,Z,g; q = c,u) |3|, Ij], 



can be significant 
mode 

.0 



y enhanced in some new physics models like the minimal supersymmetric 



(MSSM) J5|, the Top color-assisted Technicolor Model (TC2) |6| and other models (cf. 
Ref.[7|). Their results make the observation of any FCNC top quark processes a robust 
evidence for new physics. 

On the other hand, a better understanding of these channels within the SM is helpful 
for judging whether there is really new physics, i.e. to deduct the SM background from 
the experimental data at a higher confidence level such that to determine the right ranges 
for the new physics parameters. Following the top quark dominant decay channel, t — )■ 
bW~^, it has been pointed out that sizable B~ mesons can be produced via the channel, 
t — )■ |(6c)[?7,]) +cPF+ [Sl, |9|, where [n] stands for the dominant ffec]- quarkonium states via the 
velocity scaling rule of the non-relativistic QCD (NRQCD) 10|. It is interesting to show 
whether the rare SM FCNC channel t — ?■ |(c(5)[?t-]) + QZ^ can also be a useful supplement 
for heavy quarkonium production, where Q stands for the heavy quark c or b. As a short 
notation, we represent the two heavy quarkoniums, i.e. the (cc)-charmonium and (cb)- 
quarkonium, as (c(5)-quarkonium. 

In the NRQCD framework, a doubly heavy meson is considered as an expansion of 
various Fock states. In addition to the production of the two color-singlet S-wave states 



\icQ) 

rule 



^S'o]i) and | (cQ) [^Si]i), a naive order estimation basing on the NRQCD velocity scaling 



10| shows that the production of the four color-singlet P-wave states |(c(5)[^-Pi]i) and 



(c(5)P-Pj]i) (with J = 0,1,2) together with the two color-octet components |(c(5)["'^5'o]8( 



and \{cQ)[^ Si]sg) shall also give sizable contributions. Here the thickened subscripts of (cQ) 
denote the color indices, 1 for color singlet and 8 for color-octet; the relevant angular mo- 
mentum quantum numbers are shown in the parentheses. It would be interesting to study 
various Fock states' contributions to make a sound estimation on the production of heavy 
quarkonium through the channel t — )■ |(c(Q)[n]) + Z^Q, and hence to be a useful reference 
for experimental studies. 

The production of heavy quarkonium itself is very useful for testing perturbative 
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QCD |ll|-|l4j. For example, in addition to the 'direct' hadronic production of Be meson 
19| . it has been shown that its 'indirect' production channels via the top-quark, the Z^'-boson 
and the ly^-boson decays are important supplement and enough large events can also be 
produced at the LHC |8|, |9|, |20|-|26|. To make a systematic study on the Be meson indirect 
production through the top quark decays via the rare FCNC channel is one of the purpose 
of the present paper. 

The organization of the paper is as follows: In Sec. II, we describe our calculation technol- 
ogy for the top quark semi-exclusive decays to the heavy quarkonium via FCNC. Then we 
present numerical results and make some discussions on the properties of the heavy quarko- 
nium production through the top quark decays in Sec. III. The final section is reserved for a 
summary. 

II. TOP QUARK DECAYS VIA FCNC TO THE HEAVY QUARKONIUM 



According to the NRQCD factorization formula 10[, the decay width for the heavy 
quarkonium production via FCNC, i.e. t{k) — > \{cQ)[n]){q3) + Z^{q2) + Q{qi), can be 
factorized as 

dV = J2dT{t -^ {cQ)[n] +Q + Z°)(0^(n)), (1) 

n 

where the non-perturbative matrix element {0^{n)) describes the hadronization of a per- 
turbative (cQ) pair with quantum number [n] into the observable hadronic state H. The 
color-octet matrix elements are to be determined experimentally, which are smaller than the 
color-singlet matrix elements by certain v"^ order and can be estimated by the NRQCD ve- 
locity scaling rules. As for the color-singlet components, the matrix elements can be directly 
related to the wave functions at the origin for the S'-wave states or the first derivative of 
the wave functions at the origin for the P-wave states, which can be computed through the 



potential models, cf. Refs.|27H32|. 
The short-distance decay width 



dtit ^ icQ)[n] +Q + Z')) = _^|Mpci$3, 



(2) 



where J2 means that we need to average over the spin and color states of the initial particles 
and to sum over the color and spin of all the final particles. 

In the top quark rest frame, the three-particle phase space can be written as 






The 1 — 7- 3 phase space with the final massive particles can be found in appendix of Refs. [9|, 
21| . Here the parameters k and qf are four-momenta of the corresponding particles. 



I(CQ)W>(?3) 



\{cQ)[n]Kq,) 



0(?32) 




Q{qi) 



FIG. 1: Feynman diagrams for t{k) — )• \{cQ)[n]){q-i) + Z^{q2) + Q{qi), where |(cQ)[n]) stands for the 
eight quarkonium Fock states \{cQ)[^ So]i) , |(cQ)p5i]i), \{cQ)[^ Pi]i) , \{cQ)[^ Pj]i) , \{cQ)[^ So^g) 
and \{cQ)['^ Si]sg) , respectively. 

As a short notation, we take t{k) — )■ \{cQ)[n]){q3) + Z^{q2) +Q{qi) to stand for the semi- 
inclusive channels: t — )■ |(cc)[n]) -|- cZ^ and t -^ \{ch)[n\) + hZ^ . The decay-width of these 
processes can be written as: 



dV 



2567r3rra? 



Y,\M\^ds^ds2, 



(4) 



where m^ is the top quark mass, and the invariant masses Si = (gi -|- q^Y ^"^^ ^2 = (?i + Qzll; 
The integration over si and S2 can be done with the help of the VEGAS program (33 1 
whose improved version can be found in the programs BCVEGPY 19| and GENXICC [34 1- 
The symbol J2 means that we need to average over the color and spin states of the initial 
particles and to sum over the color and spin of all the final particles. With the help of the 



formulas listed in |9|, |2l|, one can not only derive the whole decay width but also obtain 
the corresponding differential decay widths that are helpful for experimental studies, such 
as dV /dsi, dT/ds2, dT /dcosOis and dT/dcos9i2, where Si = (gi + q^)^, S2 = {qi + ^2)^, ^'13 
is the angle between gl and q^, and 6*12 is the angle between gi and g*2- 

Our remaining task is to deal with the hard-scattering amplitude for the specified pro- 
cesses 

t -^ {cc)[n] + cZ° and t -^ {cb)[n] + hZ^ . (5) 

These amplitudes can be generally expressed as 

m 

iMss' = C Usi{qi) Y^ AnUs'jik), (6) 

n=l 

where m = 2 stands for the number of Feynman diagrams of these two processes, s and s' are 
spin states, i and j are color indices for the outing Q quark and the initial top quark, respec- 
tively. The Feynman diagram of the process is presented in Fig.([T]), where the intermediate 
gluon should be hard enough to produce a cc pair or a bb pair, so the amplitude is pQCD 
calculable. The overall factor C = C^ or Co stands for the specified quarkonium in color- 
singlet and color-octet, respectively. C, = .,^'£^g^^ S,j and Co = ^Jfl^^ {V2T^T''T%j, 
where ^J2T^ stands for the color of the color-octet quarkonium state and Qw is the Wein- 
berg angle. The amplitude An for each hadronic state can be read out from these Feynman 
diagrams and we put them in the appendix. 

Because of the emergence of the massive-fermion lines, the analytical expression for the 
squared amplitude becomes too complex and lengthy. We adopt the 'improved trace tech- 
nology' to simplify the amplitudes M^gt at the amplitude level. In different to the helicity 
amplitude approach, which is to get the numerical (complex) value for the whole amplitude 



at the amplitude level [19 



35I437J] , only the coefficients of the basic Lorentz structures are 
numerical at the amplitude level for the 'improved trace technology'. Because the basic 
Lorentz structures are limited, by using the 'improved trace technology', one can get the 
analytic expressions for the amplitude and conveniently result in the numerical value for the 
squared amplitude, thus, the numerical efficiency can also be greatly improved. 



24|. 



The 'improved trace technology' has been suggested and developed by Refs. [9|, I20I - 
As an explanation of the approach, we first arrange the amplitude Mss' into four orthogonal 
sub-amplitudes Mi according to the spins of the outgoing quark Q with spin s and the 
top quark with spin s', then transform these sub- amplitudes into a trace form by properly 



dealing with the massive spinors with the help of an arbitrary light-like momentum ko and 
an arbitrary space-like momentum ki, which satisfies kf = —1 and ko ■ ki = 0. The final 
results should be independent of ko and ki, which provides a way to check the rightness 
of the derived results. One can choose them to be those that can maximally simply the 
amplitude. Then we do the trace of the Dirac 7-matrix strings at the amplitude level, which 
finally result in analytic series over some independent Lorentz-structures. 

Following the standard procedures of the 'improved trace technology', we can derive 
the independent Lorentz-structures together with their coefficients for each hadronic state 
|(c(5)['T']). Here to short the paper, we will not present them here. The interesting reader 
may consult Refs. |9|, |20|-|24|| for detailed technology. The mathematical program for these 
production channels are available upon request. 



III. NUMERICAL RESULTS 

We adopt the following values to do the numerical calculation: rriz = 91.1876GeV, 
rrit = 172.0GeV, nic = l.SOGeV and nib = 4.90GeV. Leading order «« running is adopted 
and we set the renormalization scale to be 2mc, which leads to as{2mc) = 0.26. The 
color-singlet non-perturbative matrix elements can be related to the wave function at the 
origin \E's'(0) = Jl/AirRsiO) and the first derivative of the wave function at the origin 



^'p(O) = ^3/AttR'p{0), and we adopt [3^ 

\Rs{cc){0)\^ = 0.810 GeV^ \R'p{cc){0)\^ = 0.075 GeV^ 

\Rs{cb){0)\^ = L642 GeV^ \R'p{cb){0)\'^ = 0.201 GeVl 

For the color-octet S-wave matrix elements, based on the velocity scaling rule and under 



the vacuum- saturation approximation, we have 



|lO 
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{{cQ)['So]s\OsCSo)\{cQ)['So]s) ^ ^s{vf{{cQ)['So]l\O^CSo)\{cQ)['So]l), (7) 



and 



{{cQ)['S^]s\OsCS^)\{cQ)[''S^]s) ^ As{v)'{{cQ)['S^]^\0,{'S^)\{cQ)['S,],) (8) 



with As{v) ~ f^ 



A. Basic numerical results 



As a reference, we calculate the decay width for the process t — )■ c + Z", which can be 

written as 

^ ^ GFm^^\p\sfic^{ew) jg 
36\/27rmtm| 



3m| Jm,2 + |p]2 + 2|p12 ^^^2 + \^2 + /^| + [^^ 



8[2cos(2^iy) + l]sin2(^VK) 



-SmcTTi^ + Snizyrnc'^ + |p|^+ 



2|pT uA^^TT^ + v/^^T^ 



(9) 

where p stands for the relative momentum between the final two particles in the top-quark 
rest frame 



mt 



\P\ 



Then, we obtain Tt^^+zo = 0.395 GeV. 



{nic - nizY] [ml - {rric + nizY] 



2mt 



t^l(cc)N) 


r (KeV) 


rt->|(cc)[n]> 
^t-fc+zO 


t -^ r]c 


70.2 


1.78 X 10"^ 


t^ J/ip 


72.0 


1.82 X 10"^ 


t ^ he 


9.01 


2.28 X 10"^ 


t -^ XcO 


7.18 


1.82 X 10"^ 


t -^ Xcl 


8.05 


2.04 X 10"^ 


t -^ Xc2 


3.00 


7.59 X 10-6 


t^\{cc)['So]s9) 


8.78u^ 


2.22^"^ X 10-^ 


t^\{cc)[^S^]sg) 


9.001;^ 


2.28?;'^ X 10"^ 



TABLE I: Decay widths and branching fractions for the charmonium production through the 
channel t — ^ |(cc)[n]) + cZ°. 

Total decay widths and their branching fractions for the channels t — )■ |(c(5)[^]) + QZ^ 
are listed in Tables [T] and [III The squared amplitude for the color-octet decay width is 
suppressed by eight times to that of the color-singlet case. As a combined effect of such color 
suppression and the relative velocity suppression (^^-suppression at the squared amplitude 
level), the color-octet channels are quite small in comparison to their corresponding color- 
singlet production channels. So in the following discussion, if not specially stated, we shall 
not include the color-octet states' contributions into our discussions. 



t^\{cb)[n]) 


r (KeV) 


rt->|(c6)[nl> 
^t^c + zO 


t^ Be 


3.54 


8.96 X 10-*^ 


t^ B* 


2.97 


7.52 X 10^*^ 


t ^ hi) 


0.20 


5.10 X 10-'^ 


t^ Xbo 


0.34 


8.71 X lO^'^ 


i-> Xbi 


0.22 


5.65 X 10^^ 


i-> Xb2 


0.02 


3.84 X 10~^ 


t^\{cb)['So]s9) 


0.44^;'' 


1.12w'^ X 10"^ 


i^|(c6)p5i]89) 


0.37u^ 


0.94w'^ X 10"^ 



TABLE II: Decay widths and branching fractions for the (c&)-quarkonium production through the 
channel t -^ \{cb)[n]) + bZ^. 
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FIG. 2: Differential decay widths dV/dsi (Left) and dr/ds2 (Right) for t -^ \{cc)[n]) + cZ°, 
where the diamond, the solid, the dashed, the crossed, the dash-dot and the dotted lines are for 

|(cc)[i5o]i), |(cc)p5i]i), |(cc)[iPi]i), |(cc)[3Po]i), |(cc)pPi]i) and \{cc)[^P2]i), respectively. 

For the charmonium production channel t -^ \{cc)[n]) + cZ^, the total decay width for 
all the P-wave states is ~ 27.2 KeV, which is about 39% (38%) of that of rjc (J/'^P)- For 
the (c6)-quarkoniuni production channel t — )■ |(c6)[r7,]) + bZ^, the total decay width for all 
P-wave states is ~ 0.78 KeV, which is about 22% (26%) of that of B^ (B*). This shows 
that the color-singlet P-wave states can have sizable contributions to the heavy quarkonium 
production. 
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FIG. 3: Differential decay width dV/dsi (Left) and dT/ds2 (Right) for t — )■ | (c6) [n]) +bZ^, where the 
diamond, the solid, the dashed, the crossed, the dash-dot and the dotted lines are for \{cb)[^ So]i) , 
|(c6)pSi]i), |(c6)[iPi]i), |(c6)pPo]i), |(c6)pPi]i) and |(c6)pP2]i), respectively. 
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FIG. 4: Differential decay width dr/dcos6li3 (Left) and dr/dcos6li2 (Right) for t -^ |(cc)[n]) + cZ°, 
where the diamond, the solid, the dashed, the crossed, the dash-dot and the dotted lines are for 

\icc)[^So]i), |(cc)p5i]i), |(cc)[iPi]i), |(cc)pPo]i>, |(cc)pPi]i) and |(cc)pP2]i), respectively. 



Further more, to show the relative importance among different Fock states, we present the 
differential distributions dT/dsi and dV /ds2 in Figs. fl2|3|) . and the differential distributions 
(iF/(icos^i3 and dT /dcos6i2 in Figs.( l4f5|) . Two invariant variables, Si = (gi + ga)^ and 
S2 = {qi + (I2Y, and 613 stands for the angle between gl and gs, 6^12 stands for the angle 
between gl and g^. The curves for \{cQ)[^ So]r) , \icQ)fSi]^), \icQ)[^Pi]i) and \icQ)[^Pj]i) 
are presented. The difference between the color-singlet S'-wave states and the color-octet 
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FIG. 5: Differential decay width dT/dcos On (Left) and dV/dcos 6*12 (Right) for t -> |(c6)[n]) + feZ°, 
where the diamond, the solid, the dashed, the crossed, the dash-dot and the dotted lines are for 
\{cb)['So]i), |(c6)p5i]i), |(c6)[iPi]i), |(c6)pPo]i), \icb)[^Pi]i) and \{cb)[^ P2]i) , respectively. 



S'-wave states is an overall color factor times the difference between the color-singlet and 
color-octet matrix elements, the shape of their curves are the same, so the curves of the 
color-octet ones are not presented in those figures. 

As for the decay channels t — )■ |(cc)[n]) + cZ^ and t — )■ |(c6)[n]) + bZ^, because the 
heavy quarks and the heavy quarkoniums are much lighter than the Z^ boson, the largest 
dr/d cos 6'i3 is achieved when the heavy quarkonium and the outgoing heavy quark moving 
in the same direction (6*13 = 0°) or the Z^ boson and the outgoing heavy quark moving 
back to back {612 = 180°), which are shown in Figs.dH [5]). This shows that the maximum 
differential decay width is obtained when the quarkonium and the outgoing heavy quark 
moving in the same direction or the heavy quarkonium and the outgoing Z^ boson moving 
back to back in the rest frame of top quark. 

Considering the LHC runs at the center-of-mass energy y/S = 14 TeV with the luminosity 
10^^cm~^s~^, we can estimate the heavy-quarkonium events through the top quark decays 
via FCNC, i.e. about 1.8 x 10"^ rjc, 1.9 x 10^ J/\l/ and 0.7 x 10^ P-wave charmonium 
events per year can be generated; about 8.3 x 10^ Be, 7.5 x 10^ B* and 2.0 x 10^ P-wave 
(c6)- quarkonium events per year can be generated. It might be possible to find J/tp and 
Be through top quark decays, since one may identify these particles through their cascade 
decay channels, J/tp — )■ fj^'^fi~ and B^ -^ J/ip + tt or P^ — )■ J/tp + ev^ with clear signal, as 
well as Z^ decays to a fair of leptons, Z^ — t- /^/^, where / stands for leptons, i.e. e, /i, r. 
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Baring the situation pointed out here and the possible upgrade for the LHC (SLHC, DLHC 
and etc. 39|) in mind, the possibihty to study the charmonium and the (c6)-quarkoniuni 
via top quark decays is worth seriously thinking about. 



B. Uncertainty analysis 

In this subsection, we shall discuss the uncertainties for the charmonium and the (cb)- 
quarkonium production through the top quark decays. At the present leading-order cal- 
culation, their main uncertainty sources include the non-perturbative bound state matrix 
elements, the renormalization scale fiR and the quark masses rrit, rrib and rric. The scale 
uncertainty can be suppressed by a higher-order perturbative calculation, or be solved by 
a proper scaling setting methods such as the recently suggested principle of maximum con- 
formality (PMC) 40|, in which the most important thing is to pick out the non-conformal 
/3-series from higher order terms and to absorb them into the effective PMC scales, thus, 
the resulting series is conformal which is renormalization scheme and scale independent. 

In the following, we shall concentrate our attention on the uncertainties caused by rrit, 
rrtb and ttIc. We let them within the range of rrtt = 172.0 ± 4.0 GeV, rrih = 4.90 ± 0.40 GeV 
and rric = 1.50 ± 0.25 GeV. When discussing the uncertainty caused by one parameter, the 
other parameters are fixed to be their central values. 



mc(GeV) 


1.25 


1.50 


1.75 


r|(cc)[iSo]i>(KeV) 


122.2 


70.2 


43.9 


r|(cc)[3Si]i>(KeV) 


125.5 


72.0 


44.9 


r|(cc)[P-u;ai>e]i>(KeV) 


68.2 


27.2 


12.5 


r|(c5)[i5o]i>(KeV) 


3.54 


3.54 


3.53 


r|(c5)[3Si]i>(KeV) 


2.92 


2.97 


3.01 


'^\{cb)[P-wave]i)(^^'^) 


1.05 


0.78 


0.62 



TABLE III: Uncertainties for the decay width of the process t — t- |(cc)[n]) + cZ^ and t — t- |(c6)[n]) + 
bZ^ by varying rric G [1.25, 1.75] GeV, where |(cc)[P — waveli) and |(c6)[P — wave]i) stands for 
the sum of the four color-singlet P-wave states for the (cc) and (cb) quarkoniums respectively. 

Typical uncertainties for rric, f^b and rrit are presented in Tables UTTl II VI and IV l respectively. 
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mfe (GeV) 


4.50 


4.90 


5.30 


r|(c5)[i5o]i>(KeV) 


4.65 


3.54 


2.74 


r|(cS)[35iU>(KeV) 


3.95 


2.97 


2.28 


r|(c6)[P-«,ai>e]i>(KeV) 


1.06 


0.78 


0.59 



TABLE IV: Uncertainties for the decay width of the process t — )• (c6)[n] by varying nib G [4.50, 5.30] 
GeV, where |(c6)i[P — wave]) stands for the sum of the four color-singlet P-wave states. 



mi (GeV) 


168.0 


172.0 


176.0 


r|(cc)[iSo]i>(KeV) 


63.9 


70.2 


76.8 


r|(c5)p5i]i)(KeV) 


65.6 


72.0 


78.7 


r|(cc)[P-u;at;e]i>(KeV) 


25.3 


27.2 


29.2 


r|(e5)[i5o].)(KeV) 


3.18 


3.54 


3.91 


r|(c6)[35i]i>(KeV) 


2.68 


2.97 


3.27 


'^\{cb)[P-wave]i)0^^^) 


0.72 


0.78 


0.84 



TABLE V: Uncertainties for the decay width of the process t — )• (cc)[n] and t — t- {cb)[n] by varying 
mt £ [168.0, 176.0] GeV, where \{cc)[P — waveli) and |(c6)[P — wave]i) stands for the sum of the 
four color-singlet P-wave states for (cc) and (cb) quarkonium accordingly. 



It shows that sizable uncertainties can be found for varying rric, mi, and rrit, respectively. 
The decay width will decrease with the increment of rric, mi, and mj, and such tendency 
slowly down with a heavier quark mass. Moreover, one may observe that the decay widths 
for P-wave states are more sensitive to the quark masses than the case of S-wave states. 
Adding all the uncertainties in quadrature, we obtain 



ri(cc)[i5o]i> 
r|(cc)[35i]i> 

'- |(cc)[P— -iiiatieji) 
\{cc)[S —'wave\sg) 



70.2+1;° KeV, 



-^20+53-6 



-27, 



?KeV, 



27.2+^°:? KeV, 



(10) 



for t — )■ |(cc)[r7,]) + cZ^] and 



l(c5)F5o]i> 



3.54+^-;^^ KeV, 



-2.38 
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|{cb)p5i]i> 



\{cb)[P—wave]i) 



\(cb)[S—wave]sg) 



2.97+|95 KeV, 
0.78+J;^^ KeV, 
0.81+°:^6 X v^ KeV. 



(11) 



Fort^ \{cb)[n]) + bZ'^. 

If assuming the high excited heavy-quarkonium states decay to the ground color-singlet 
and spin-singlet S wave state |(c(5)[^>S'o]i) with 100% efficiency via the electromagnetic or 
the hadronic interactions, then we obtain the total decay width of top quark decay channels 



(t^|(cc)[iSo]i>+cZO) 



(t^|(c5)[i5o]i>+6Z") 



171.1+^f/KeV 
7.32+^-49 KeV 



(12) 
(13) 



where we have set v'^ = 0.3 for the color-octet S-wave charmonium states |(cc)[^S'o]8fi') and 
\{cc)[^Si]sg), and v"^ = 0.2 for the color-octet S-wave (c6)-quarkonium states \{cb)[^ So]sg) 
and \{cb)[^Si]8g), respectively. 
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FIG. 6: Uncertainties of dV/dsi (Left) and dT/ds2 (Right) for t — t- |(cc)[n]) + cZ^, where the 
contributions from the color-singlet S-wave and P-wave states have been summed up. The solid, 
the dashed and the dotted lines are for nic = 1.25 GeV, 1.50 GeV and 1.75 GeV, respectively. 

To show how the decay widths depend on the quark masses, taking t — )■ |(cc)[?t,]) -I- cZ° 
as an example, we present its differential decay width with several typical rric and mb in 
Figs. ipiTl) . Uncertainties of dV/dsi and dT/ds2 of the process are presented in Fig.(l6l), 
where the contributions from the color-singlet S'-wave and P-wave states have been summed 
up; Uncertainties of (iF/ci cos 6'i3 and ciF/ii cos 6*12 of the process are presented in Fig. ([7]). In 
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FIG. 7: Uncertainties oi dT/dcosOi^ (Left) and dT/dcos6i2 (Right) for t — )• |(cc)[n]) +cZ^, where 
contributions from the color-singlet S'-wave and P-wave states have been summed up. The solid, 
the dashed and the dotted lines are for rric = 1.25 GeV, 1.50 GeV and 1.75 GeV, respectively. 

these figures, the contributions from the color-singlet S'-wave and P-wave states have been 
summed up. 

IV. CONCLUSIONS 

Being by far the heaviest fundamental particle and the only quark decaying before 
hadronization takes place, the top quark provides a privileged tool for precise studies of 
the SM. A precise study of the top quark is helpful for understanding the hadronic and 
electroweak interactions and for searching new physics beyond the SM. 

We have made a detailed study on the heavy-quarkonium production through the top 
quark semi-exclusive decays via FCNC, t — )■ \{cQ)[n]) + QZ^, within the NRQCD framework. 
Results for the eight quarkonium Fock states, i.e. |(c(Q)[^S'o]i,8), |(c(5)P'S'i]i,8), |(cO)[^-Pi]i) 
and |(c(5)[^-Pj]i) have been presented. To provide the analytical expressions as simplify as 
possible and to improve the numerical efficiency, we have adopted the 'New Trace Technol- 
ogy' for dealing with the hard scattering amplitude directly at the amplitude level. 

Numerical results show that the color-singlet P-wave states in addition to the S-wave 
states can also provide sizable contributions to heavy quarkonium production, so one also 
needs to take the P-wave states into consideration for a sound estimation. For the charmo- 
nium production channel t — )■ |(cc)[n]) +cZ^, the total decay width for all the P-wave states 
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is ~ 27.2 KeV, which is about 39% (38%) of that of rjc (J/'ip)- For the (c&)-quarkonium 
production channel t — )■ |(c6)[r7,]) + bZ^, the total decay width for all P-wave states is 
~ 0.78 KeV, which is about 22% (26%) of that of B^ (B*). If assuming the high excited 
heavy-quarkonium states decay to the ground color-singlet and spin-singlet S wave state 
|(c(5)[^5'o]i) with 100% efficiency via the electromagnetic or the hadronic interactions, then 
we obtain the total decay width of top quark decay channels via FCNC 



t^|(c5)[i5o]i> 



171.1+^f/KeV, 



7.32+^-49 KeV. 



ii^|{c6)[iSo]i> 

At the LHC, due to its high collision energy and high luminosity, sizable heavy-quarkonium 
events can be produced through top quark decays via FCNC, i.e. about 10^ rjc and 10'^ Be 
events per year can be obtained. 
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Appendix A: An for t -^ {cQ)[n] + QZ^ 

For convenience, we present the An for the hard scattering of the processes t{k) — )■ 
\{cQ)[n]){q3) + Z^{q2) + Q{qi) based on the Feynman diagrams shown in Fig.(Il]), where Q 
stands for c or 6 quark respectively. 

For the (c(5)-quarkonium in S-wscve states An can be written as 



A2 



■ ny)(g) 



rrir 



7/. 






m- 



-,iiq2){Tc-ir>) 



Jg=0 



(g32 + gi)^ 



/(g2)(Tc-75) 



'31 



rrit 



(^2 + gsi) 



m\ 



r7M 



9=0 



For the (cQ)-quarkonium in P-wave states An can be written as 



(Al) 
(A2) 



^? 



S=0,L=1 



e^q^ 



dqa 



7/. 



Kii) ., ^1 + ^, 



TTlr 



7m^ 



(g32 + gi)^ (gi + 53) 



mt 



(g2)(Tc-75) 



9=0 



(A3) 
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A. 



S'=0,L=1 



^r(?3 



dqa 



■'.s^jMiT.-..) /j:s:'^:. -r. 



'(g32 + gi)^ 



mi 



At 



S=1,L=1 



^a/3(?3)3- 



A. 



S=1,L=1 



^iMi) 



d 



7m 
7m 



<(^) .^ ^1+^3 + "^ 



77^: 



KM) 



{(12 + gsi) 

/(g2)(Tc-75) 



g=0 



Jg=0 



(fe + gi) 



.)(^-7a) /'^V":. ^. 



(g2 + qsi] 



rat 



(A4) 
(A5) 
(A6) 



q=0 



Here T^ 



|sin^(6'vi/). e(g2) is the polarization vector of Z^ boson, ^^(gs) and 5;(g3 



are the polarization vectors relating to the spin and the orbit angular momentum of {cQ)- 
quarkonium, s^^Jq^) is the polarization tensor for the spin triplet P-wave states with J = 0, 1 
and 2, respectively, gsi and q^2 are the momenta of the two constitute quarks, i.e. 



gsi 



m. 



-gs + g, g32 



mq 



Q3-Q- 



(A7) 



where q stands for the relative momentum between the two constitute quarks in (cQ)- 
quarkonium. M ^ rric+mQ are adopted to ensure the gauge invariance of the hard scattering 
amplitude. The projectors n°g(g) and n^g(g) are for spin-singlet and spin-triplet states 
accordingly, and their covariant form can be written as 41 1 

-\/M 



KM 



ArricmQ ^^32 

-Vm 



(^32-"^q)75(^31 + "^c) 



AnicmQ '^32 
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